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We calculate the nonequilibrium mean-field 'temperature' of a Brownian system 
in contact with a heat bath. We consider two different cases: an equilibrium bath 
in the presence of strong external forces and a nonequilibrium bath. By proving the 
existence of a generalized fluctuation-dissipation relation this mean-field 'tempera- 
ture' can be used to describe a nonequilibrium system as is if it were in thermal 
equilibrium with a thermal bath at the mean-field 'temperature' mentioned above. 
We apply our results to chemical reactions in the presence of external forces showing 
how chemical equilibrium and Kramers rate constants are modified by the presence 
of these forces. 
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I. INTRODUCTION 

In the thermodynamic analysis of mesoscopic and macroscopic systems out of equilibrium 
an interesting question arises: Can we define a nonequilibrium 'temperature' ? Although this 
is possible, this 'temperature' cannot be a thermodynamic temperature since there is not a 
thermodynamic zero principle behind it. The consequences of this fact have been studied in 
detail in Ref . p] . Nevertheless, the concept of a nonequilibrium 'temperature' can be used to 
parametrize the quasi-equilibrium states of a nonequilibrium system PQ , [2] , [3] , [Ij , • 

The notion of a nonequilibrium 'temperature' can be made more precise by noting that 
this is a statistical concept related to the energy of a system involved in the erratic motion 
of its Brownian degrees of freedom, that is, related to the thermal energy of the system. 
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Since the thermal energy of equihbrium and nonequihbrium systems will in general differ, 
one expects that the nonequihbrium 'temperature' of an out of equilibrium system will be 
different from that of the heat bath. Introducing the nonequihbrium 'temperature' has the 
advantage that one may describe the system as is if it were at thermal equilibrium with a 
hypothetical bath with a temperature corresponding to this nonequihbrium 'temperature'. 
This possibility has been studied for example for quantum corrections to the low temperature 
of thermodynamic systems or thermal radiation not in equilibrium in Ref. |6] (page 104, 
section 34 and page 189, section 63, respectively) and in the case of granular matter in 
Ref. Here, by using the nonequihbrium 'temperature', we will show that in the case of 
chemical reactions an external force shifts chemical equilibrium and increases the reaction 
rates. This effect may be particularly important in the case of photochemistry where the 
incident light increases the velocity of the reaction. 

In nonequihbrium systems such as small, confined or glass- like systems jTj, introducing 
a nonequihbrium 'temperature' is appropriate since this concept is a consequence of the 
existence of internal constraints or long-range forces and correlations which maintain the 
system out of thermodynamic equilibrium. Some of these constraints may appear at low 
temperatures and for small masses as in the case of quantum effects [6j mentioned above. In 
other cases, it is a purely classical effect when the range of the interactions is similar to the 
size of the system. 

In this article, we propose a definition of the nonequihbrium 'temperature' T(a, t) in 
analogy with the equipartition theorem 



where ks is Boltzmann's constant, f{b, a,t) is a probability distribution describing the state 
of the system which, for convenience's sake, we assume to depend on a pair of conjugated slow 



ity distribution. One of the goals of this article consists of showing how this nonequihbrium 
'temperature', which we will call effective 'temperature', is connected to the bath tempera- 
ture and the forces which maintain the system out of equilibrium |3]. We also analyze how 
the fluctuation-dissipation theorem (FDT) is modified in systems far from equilibrium and 
how the effective 'temperature' comes into play in order to extend the validity of this FDT 
to quasistationary nonequihbrium systems. 





a and fast b variables (a velocity). In addition, p{a,t) = J f{b,a,t)db is a reduced probabil 
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The paper is organized as follows: In section 2 we analyze Brownian motion in a field 
of force. We obtain the Fokker-Planck equation and a general expression for the effective 
'temperature' and its average value the mean-field 'temperature'. We also obtain the gen- 
eralized Smoluchowski equation describing the quasi-equilibrium state. Section 3 is devoted 
to studying the effects of the mean-field 'temperature' on chemical reaction rates. In section 
4 we analyze the FDT. Finally, in section 5 we present our main conclusions. 



II. BROWNIAN MOTION IN A FIELD OF FORCE 

Let us consider a one dimensional Brownian gas in contact with a heat bath at temperature 
To [8]. Let be H{r) = mv? /2+G{x) the Hamiltonian of a particle of mass m, where F = (x, u) 
represents a point in the one-particle phase space and G{x) being the external potential. 

The analysis of this Brownian gas is based on thermodynamics through the definition of 
entropy by means of the Gibbs entropy postulate P| 

S{t) = -kB [ /ln-^rfF + 5eq., (2) 

J Icq. 

where ks is Boltzmann's constant, f{T,t) is the phase-space distribution function, Seq. is 
the equilibrium entropy and 

/eq.~exp(-^l, (3) 



is the equilibrium distribution function. Variations in the probability density f(T,t) cause 
changes in the entropy which can be obtained from Eq. ^ 

6S = - [ (ksln^ + ^^)6fdr, (4) 



fcq. Tq 

where we have taken into account that SSeq. = — J (/igq./To) 6fdT. Here /ieq.(^5 1) = —G{x) + 
fiQ is the equilibrium chemical potential (mechanochemical potential) per unit of mass and 
^q{p,Tq) is the corresponding thermodynamic potential, with p being the pressure. By 
defining the nonequilibrium chemical potential 

/ 

/i(F,t) = ksTohl- h/ieq., (5) 

Jeq. 

the thermodynamic quantity conjugated to the density f{T,t), it is possible to write Eq.(|4]) 
in the compact way 

To6S=- J fi{T,t)5fiT,t)dT (6) 
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which constitutes the Gibbs equation in phase space. 

A gradient of the chemical potential in phase space ([s]) induces a diffusion process which 
tends to restore the equilibrium state. Through this process, the distribution function 
changes according to the generalized Liouville equation 

-fir, t) + v{T, t) ■ Vr/(r, t) = -^^(r, t) (7) 

which defines the diffusion current J{T,t) and where v{T,t) = {x,u) = {u, —VG/m) is the 
velocity corresponding to the hamiltonian flow and Vr = (V, d/du), with V = d/dx. 
From Eq. ^ we can obtain the rate of change of the nonequilibrium entropy 



which in combination with Eq. ([T]) and after partial integration leads to 

f = ^ {{-VG)u) - ^ I j(r,t)|-/i(r,t)rfr. (9) 

Here, the first term on the right hand side of Eq. ^ is given by 

{i-VG)u) = j [-VGix)]ufiT,t)dT (10) 

and constitutes the power supplied by the external field which is dissipated in the system 
as heat. This contribution to the entropy change comes from the Hamiltonian evolution 
of the distribution function and therefore cannot be assimilated into the contribution due 



to diffusion. Thus, Eq. (10) can be interpreted as the rate of heat exchanged with the 
surroundings 

f^ = {{~VG)u), (11) 

with dq being the amount of heat released in a time dt. In addition, the second term on 
the right-hand side of Eq. ^ constitutes the entropy production rate due to irreversible 
processes a which must be non-negative (a > 0) according to the second law. Therefore, the 
rate of change of entropy can be written in the compact way 

dS _ ^ dq ^ ^ ^^2^1 
dt To dt ' 

expressing the balance between the exchange of heat with the surroundings and the entropy 
generated in the irreversible processes established in the system. The entropy production 
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contains the current J{T,t) and its conjugated thermodynamic force {d / du) ^{T , t) . Fol- 
lowing the postulates of nonequilibrium thermodynamics [9], these quantities are related 
through the phenomenological law 

L d 



J{T,t) 



Todu 



/i(r,t), 



(13) 



where L is the phenomenological coefficient. By using the expression of the nonequilibrium 
chemical potential, Eq. ^ and Eq. (13) one obtains 



+ u\fiT,t), 



(14) 



m du 

where we have identified L/ J'Tq as the friction coefficient ( of the Brownian particle {L/J'Tq = 
C/m). Hence, from the definition of a given through Eq. ([o]), along with Eq. (14) we obtain 

mC f [{kBTo/m) (d/du) /(E, t) + /(F, t)u]' 
To 



a 



-dr. 



(15) 



/(r,t) 

It is worth to emphasize that a stationary nonequilibrium state {dS/dt = 0) exist for which 

To^ = ^''^ 
This corresponds to a stationary state of nonzero entropy production which differs from the 
equilibrium state characterized by die^q/dt = aieq = 0, a condition which is satisfied by the 
local Maxwellian (|3|. 

Finally, by substituting Eq. (14) into Eq. ([T]) we obtain the Fokker-Planck equation 
describing the dynamics of the Brownian gas 



di 



f{r,t) = -uvf + v 



G{x) 



m 



^f + c(^-^—f + —uf\ 
du \ m dv? du 



(17) 



III. EFFECTIVE AND MEAN-FIELD 'TEMPERATURES' 



At this point, assuming that the velocity u is the fast variable (in Brownian motion inertia 
usually constitutes a short lived effect), it is appropriate to write 

/(F,t) = 0,(M,t)p(x,t), (18) 

where (pxiu^t) is the conditional probability density and p{x,t) = mj f(T,t)du is the con- 
figurational probability density which evolves according to 

d 



dt 



p{x,t) = -V 



m / uf{T,t)du 



(19) 
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Equation ( 18 ) expresses the coupling between the macroscopic process triggered by the field 



G{x) and the microscopic process of the momentum relaxation [TO]. Equation (19) has 



been obtained by partial integration of Eq. (17) over velocity and thus implicitly defines 



the diffusion current J{x,t) ^ m J uf{T^t)du = p{x,t)v{x,t). This current satisfies the 
evolution equation 



Ft 



J{x,t) + CJ{x,t) = p{x,t) 



-V 



G{x) 



m 



- V 



u^(j)x{u, t)du I p{x, t) 



(20) 



obtained from Eq. (17) by multiplying by u an integrating by parts with the appropriate 



boundary conditions. Here, the second moment of (j)ri.{u,t),J u^(j)x{u,t)du is proportional to 
the thermal energy, i. e. the amount of energy of a system necessary for the erratic motion of 
its Brownian degrees of freedom. This suggests the definition of an effective 'temperature' 



ksT {x,t) 



m 



u^<pxiu, t)du 



(21) 



in analogy with the equipartition theorem [3],[1],[5]. If in particular, the initial distribution 
is given by 

(22) 



/(r,0) = p(x,0)exp { -- (m -f(x,0)) m/kBTo 



the solution of Eq. (17) at later times will have also the same Gaussian form and is given 



by (see Ref. [2] CH. IX, §8) 



1 



/(r,t) = p{x,t)exp <{--(«- v{x,t)) m/kBTo [> . 



(23) 



In this case the second moment of (f)x{u,t) becomes 



m 



+ v{x,ty 



(24) 



which defines the effective 'temperature' 



ksT (x, t) = fc^To + mv{x, t^. 



(25) 



This effective 'temperature' enters the expression of the diffusion current which, from Eq. 



(20) for long times as compared to C \ reduces to 
J{x,t) = -D{x,t) 



Vp(x,t)+ Z^';'^^ V1>(x,t) 
kBi [X, t) 



(26) 
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where now 



G{x) + kBT{x,t) is an effective potential and D{x,t) 



{kBT{x,t)/m)( ^ is the bare effective diffusion coefficient. Substitution of Eq. (26) into 



(19) yields the generalized Smoluchowski equation 
d 



dt 



p{x,t)=V {D{x,t) 



Vp(x,t) + -4^V$(x,t) 
kbI [x, t) 



(27) 



If the conditions are such that the system is in a quasistationary state characterized by 



Eq. (16), thus, from Eqs. (11), (15) and (|23]) one obtains 

'G{x) 



v{x,t) = C"V 



m 



which implemented in Eq. ( 25 ) gives 

kBT(x,t) = fcsTo + — \C-^VG(x 
m 

A particular simple case corresponds to the potential G{x) 
reduces to 



(2^ 



(29) 



-Fx for which Eq. (29) 



kBT (t) = kBTo + -{(-' Py. (30) 

Another interesting particular case corresponds to the presence of thermal forces. The 
existence of, for example, an homogeneous temperature gradient in the bath brings the 
system to a nonequilibrium state. This non-homogeneous bath temperature Tnh(x) yields 
a Fokker-Planck equation in which G{x) = 7lnTnh with 7 a characteristic energy, [llj and 
therefore 



1 



kBT {x, t) = kBT^h + - V In r„h) 



(31) 



Eqs. (29)-(31) show that when the system is subjected to a large force, the averaged 
kinetic energy related to Brownian motion increases, enabling one to define an effective 
'temperature' differing from the corresponding temperature of the bath. 

In addition, for practical purposes it is convinient to introduce an estimate of the 'tem- 
perature' field T{x,t) which allows to incorporate the effects of the strong external forces 
in such a way that the usual techniques of nonquilibrium statistical physics can be used to 
calculate the properties of the system. In the lowest approximation, this can be done by 
means of the mean value {T{x,t)) = J T{x,t)p{x,t)dx, which we will call the mean-field 
'temperature' T^f-. 



kBT^iit) = kB (To) + -{[C'VG{x)Y). 



m 



(32) 



8 



Adopting this mean-field 'temperature', the current (26) takes the form 

J(x,t) = -D 



Vp(x,t) + ^VG(x) 



which, after substituted into Eq. (|19j), yields the generalized Smoluchowski equation 

d 



dt 



p{x, t) = D^iV^p{x, t) + C"' V t)VG{x)] 



(33) 



(34) 



where Dmf = {kBTtn{/fn)C ^ is the mean-field diffusion coefficient. At quasi-equilibrium, the 
solution of Eq. ( 34 ) 



exp 



G{x) 



knT, 



B-i mf 



(35) 



is characterized by a mean-field thermal energy ksT^^. In Eq. (35) we can interpret the 



mean-field 'temperature' as the 'temperature' for which the configuration probability density 
Pqo of the nonequilibrium system is equal to that given by Boltzmann's probability distribu- 
tion formula for a equilibrium system with 'temperature' T^i- It is convenient to notice here 



that the previous results (33), (34) and (35) are exact when a constant force is applied on 



the system, because in this case Eqs. (30) and (32) are equivalent. 



From the previous analysis, it is plausible to assume that inhomogeneities caused by the 
application of forces of different type may introduce similar effects in the dynamical prop- 
erties of the system. For example, this will be the case of diffusion in the presence of a 
shear flow [I2] , [13] , [13] • Finally, it is worth stress that introducing the mean- field 'tempera- 
ture' Tmf is consistent only if A(mf^) << ksT^. Otherwise, the complete nonhomogeneous 
nonequilibrium effective 'temperature' T{x,t) has to be considered completely. 



IV. MODIFIED CHEMICAL REACTION RATES 



As an application of practical interest of the previous analysis, let us examine how the 
application of a large external force on a chemical system modifies the conditions of chem- 
ical equilibrium. In order to perform this analysis, we will assume that x represents the 
reaction coordinate and that the reaction itself can be described as a diffusion process along 
the reaction coordinate [I5]. For chemical systems G{x) is the free energy controlling the 
reaction. 



As in Ref. [15], the current J{x,t) given through Eq. (33) can be rewritten in the form 

J{x, t) = -D^{ e-^(^)/'=^^-f Ve^-f(^'*)/'^«^-f , (36) 
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where we have introduced the chemical potential 

Atmf(a;,t) = /csTrnflnp + G'(x). (37) 

When the height of the energy barrier separating the two minima of the potential is large 
compared to thermal energy, a fast relaxation towards the local minima occurs. Then, the 
chemical potential becomes a piece-wise continuous function of the coordinates 

fJ^mfix, t) = fx^fix, t)9(xo -x)+ /i+f(x, t)Q{x - Xo), (38) 

with /i~£ and yU^f referring to the chemical potential at the left and right wells, respectively. 
Consequently the probability density also splits as 

p(x, t) = pi(t)e-[^(")'^("^)]/'=^^-^e(xo -x)+ p2(t)e~[^(")~^("^^l/'^'^-'e(x - xq), (39) 

here pi{t) = p{xi,t) with i = 1,2 are the values of the probability density at the minima, 
0(x) is the step function and Xj with j = 0, 1,2 are the coordinates of the maximum and 
the minima of the potential, respectively. 



Starting from the mean-field Smoluchowski equation (34) and using relations (36)- (39), 



it is possible to derive the following kinetic equation for the concentrations pi 

kfp2-k^'pi, (40) 



dpi dp2 f f 



= D^f^ exp 



(41) 



dt dt 

where the forward and backward reaction rates are given by 

' G{x2,i) - GixoY 
ksTtaf 

From this equation it follows an interesting result that can be expressed as the ratio between 
the mean-field rate fc™^ and the reaction rate when no external field is applied k: k™^/k = 
exp (AG/kBTo) exp[— (AG/kBTo) e] where e = To/T^f is a measure of the deviation from 
thermal equilibrium and AG = G{xo) — G{x2,i) that is a positive quantity by definition 
[16], [T7]. Figure 1 shows k^^/k as a function of e. It is clear that the external force shifts 
the chemical equilibrium and increases the reaction rates. This is particularly interesting 
in photochemical reactions because the correction involved in the mean- field 'temperature' 
is proportional to the energy of the electromagnetic field, that is to the number of photons 
associated to the incident light. Another interesting case is when the chemical reactions 
take place in a nonequilibrium medium like alive organisms or chemical reactors in which 
the presence of thermal or concentration gradients as well as hydrodynamic flows modify the 
chemical equilibrium. 
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Thermal energy: k T , (J) 



FIG. 1: Schematic representation of the mean-field reaction rate k^f as a function of the temper- 
ature of the bath for different values of the correction term a = ^{{C~^'^G)^) with AG = 1 and 
Dmf\/ G" {x2,i)\G" {xo) \ / {2TTkBTj^{) = 1. The inset shows the ratio k'^^/k of reaction rates as a 
function of e = To/T^f for three different values of AG/ksTo. When the force associated to the 
energy G{x) is large the thermal energy is large enough to increase the mean-field reaction rate. 

V. GENERALIZED FLUCTUATION-DISSIPATION THEOREM 



Once the generalized Smoluchowski equation (34) has been formulated, it is convenient 
to derive the generalized fluctuation-dissipation theorem which relates the time derivative 
of the correlation function CAF{t,t') with the response function RAp{t,t') of an observable 
A{x) through the mean-fleld 'temperature' when at time t = to the system is perturbed by 
the external fleld F{x,t) [5] , [18] , [19] . As a consequence of this and in analogy with slow 
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relaxing systems p], the generalized fluctuation-dissipation theorem offers a statistical tool 
in order to evaluate the mean-fleld 'temperature'. 

For a sufliciently weak perturbation around the quasi-equilibrium state characterized by 



Eq. (35) and the mean-fleld 'temperature' T^f, the response of the system can be expressed 



in terms of the deviation of the average value {A{t))p in the presence of the force fleld F 
with respect to the unperturbed case {Ait))^^- 

{A{t))F-{A{t))^,= [ RAF{t,t')F{t')dt', (42) 



to 



where we have deflned 



A{x)pF{x,t)dx and {A(t)) 



qe 



A{x)pqc{x, t)dx. 



(43) 



To obtain an explicit expression for the left-hand side of Eq. (43), we may use the fact 
that the perturbed probability density pp satisfles the identity 



pp{x,t) = J GF{x,t\x',t')pF{x',t')dx', (44) 
where GF{x,t\x' ,t') is the Green function in the presence of the perturbation. Considering 



also that the quasi-equilibrium solution of Eq. (34) in the presence of the perturbation F is: 



qe 



PqeZF{x,t), with Zpi^Xjt) = Z[F{x,t)] a functional of the perturbation fleld F. Now, 



after substituting p^^ in Eq. (|44|) and rearranging terms one obtains 



Pqe{x,t) = I Gf{x, t\x', t' ) ^^^^ ' ^ ^ Pqc(x , t')dx' . 

Zf[x, t) 



(45) 



Using Eq. (45) and the fact that p{x,t) = J G{x,t\x' ,t')p{x' ,t')dx' with G{x,t\x' ,t') the 
unperturbed Green function, we may establish the quasi-equilibrium fluctuation relation 



GF{x,t\x',t') = ■^^j^^^G{x,t\x',t'), 

Zjf[X , t ) 



(46) 



which expresses the perturbed Green function in terms of the unperturbed one. Now, using 
that pF{x,to) = pqe(a;,to) at t' = to, then the following relation holds 

(47) 



PF{x,t) = / GF{x,t\x',tQ)p^c{x',to)dx'. 



The substitution of Eq. (47) in Eq. (43) leads, after using the fluctuation relation (46), to 
the formula 

{A{t))F= [ f A{x)^4^;^G{x,t\x',to)p^,{x',to)dxdx'. (48) 



ZF{x',to) 
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Here, Eq. (48) may be written in a form similar to Eq. (42) by approximating 



Z[F{x,t)]/ Z[F{x' ,to)] in its series expansion up to first order in F. This operation gives the 
integral relation 



{m)F - {A{t)) 



qe 



(49) 



Recalling now that for the Smoluchowski operator 

exp{FoB{x)/kBT^i) 



Z[Fix,t)] 



(50) 



{exp{FoB{x)/kBT^f))' 

where we have used F{x, t) = FodB{x)/dx for convenience and assumed that Fq is a constant 



[8]. Therefore, using (50) in Eq. (49), for weak perturbations exp{FQB /ksT) ~ l — F^B /ksT 
and we obtain 



{A{t)), - {A{t)) 



qe 



1 



-CAB(t, to) 



(51) 



where CAB(t,to) = {A(t)B(to)) — {A(t)) {B (tQ)) . Assuming now that the quasi-stationary 



state of the system varies slowly enough, then we can replace Tmf(t) by Tmf(to) in (51) 



The resulting expression can then be written as the integral of the time derivative of the 



correlation function which, after being compared with Eq. (42) finally gives 

d 



RAsii, to) 



1 



-CAB(t, to 



(52) 



fci?Tmf(to) dto 

which constitutes the quasi-equilibrium fluctuation-dissipation relation valid for times t > to- 
This result implies that the fluctuation-dissipation relation can be generalized to the quasi- 
equilibrium state by incorporating the corrections on system's temperature due to the large 
potential V{x), |H],[2n]- Clearly, this result is fully compatible with the quasi-equilibrium 



fluctuation relation (46). 



VI. CONCLUSIONS 



In this paper we have introduced the concepts of nonequilibrium effective 'temperature' 
T{x,t) and its average value, the mean-field 'temperature' Tmf, by using a generahzation of 
the equipartition theorem. This mean-field 'temperature' parametrizes the quasi-equilibrium 
or metastable states of a nonequilibrium system. This parameter reveals that in the presence 
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of large external forces the thermal energy of the system increases, therefore promoting the 
thermal motion of its corresponding degrees of freedom. 

We have shown that at equilibrium this mean-field 'temperature' reduces to the bath tem- 
perature while at quasi-equilibrium it contains quadratic corrections coming from the forces, 
either Hamiltonian or thermal, acting on the system. The existence of internal constraints 
like surface effects or long-range forces would lead to similar phenomena as the ones arising 
from the applied external forces. 

Once this mean-field 'temperature' has been introduced, we can describe a nonequilibrium 
system as is if it were in equilibrium with a thermal bath at this T^f since it has been proven 
that fluctuations in a nonequilibrium system satisfy a generalized fluctuation-dissipation 
theorem containing T^f. 

We have illustrated the possible implications in practical situations by analyzing how 
external forces modify the chemical equilibrium and the Kramers reaction rates. We point 
out that external forces increase the reaction rates, such as in photochemical reactions or in 
the presence of nonequilibrium substrates. 
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Thermal energy: k T (J) 
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